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Abstract 

Abstract The quasi-SU(3) symmetry, as found in shell model 
calculations [0, refers to the dominance of the single particle plus 
quadrupole - quadrupole terms in the Hamiltonian used to describe 
well deformed nuclei, and to the subspace relevant in its diagonaliza- 
tion. It provides a very efficient basis truncation scheme. It is shown 
that a small number of SU(3) coupled irreps, those with the largest 
C2 values within the direct product of the proton and neutron SU(3) 
irreps with spin and 1 (for even number of particles) , and spin 1 /2 
and 3/2 for (for odd number of nucleons), are enough to describe the 
low energy spectra and B(E2) transition strengths of 21 Ne, 23 Na and 
25 Mg. A simple but realistic Hamiltonian is employed. Results com- 
pare favorably both with experimental data and with full shell model 



*Electronic address: cvargas@fis.cinvestav.mx; fellow of CONACyT. 
^Electronic address: hirsch@nuclecu.unam.mx 
■i-Electronic address: draayer@lsu.edu 



1 



calculations. Limitations and possible improvements of the schematic 
Hamiltonian are discussed. 

PACS numbers: 21.60.Fw, 21.60.Cs, 27.70.+q 

Keywords: Quasi-SU(3) symmetry, energies, B(E2) values, 21 Ne, 23 Na 
and 25 Mg. 



1 Introduction. 

Since its introduction more than fifty years ago [@, the shell model has been 
a fundamental tool in the microscopic description of nuclear properties. Full 
shell model calculations in the sd- Q and /p-shells [fl] provide very accurate 
predictions for energy levels, electromagnetic transition strengths and weak 
decay half lives. State of the art codes allow studies up to A = 50 ||. The 
system symmetries, although in general are not exact, provide a natural trun- 
cation scheme of the Hilbert space while keeping the predictive power of the 
theory, and provide a qualitative understanding of collective nuclear modes. 
Algebraic models are particularly well suited to describe systems with sym- 
metries, using either bosonic representations ||, fermionic representations 
|], 0], or combination of both ||. 

In the present article we concentrate our attention in the fermionic SU(3) 
algebraic model developed by Elliott ||. Based on the crucial role the 
quadrupole - quadrupole interaction plays in deformed systems, the SU(3) 
algebra is the natural language to describe quadrupole excitations in a har- 
monic oscillator basis in light nuclei. The pseudo SU(3) model [0], built 
over the pseudo spin symmetry, plays a similar role in heavy deformed nu- 
clei. Used as an approximate symmetry, i.e. allowing the mixing of different 
irreducible representations (irreps) through the single particle energies and 
the pairing interaction, it provides a very good description of low energy 
bands, B(E2) and B(M1) transition strengths in the rare earth region, both 
for even-even || and A-odd [|ItJ deformed nuclei. 



In the pseudo SU(3) model the truncation of the Hilbert space usually 
excludes intruder orbits, which are known to provide an important contribu- 
tion to the total quadrupole moment. The quasi SU(3) truncation scheme 
[H] offers a simple and consistent way to include intruder orbits in the SU(3) 



description of heavy deformed nuclei. In |L1]] it was shown that including 
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the leading SU(3) irreps (those with the largest Ci values) and spin and 1 
the interplay between the quadrupole - quadrupole interaction and the spin 
- orbit splitting is well described. A detailed description of four even-even 
nuclei in the s<i-shell, ranging from 20 Ne to 28 Si was presented in [pT2 ] . Three 
odd-odd nuclei in the sd-sheh 1 are discussed in The present article deals 
with the odd-mass nuclei in the s<i-shell. In the three cases the main goal is 
to prove that in well deformed nuclei, even in the presence of a large spin- 
orbit splitting, it is possible to make a good description of the low energy- 
spectra using a small number of SU(3) irreps. This is the meaning of the 
"quasi SU(3) truncation scheme". The results presented here and in the two 
accompanying articles strongly support this conclusion. They open the pos- 
sibility for a coherent description of nucleons in normal and intruder levels 
in heavy deformed nuclei using the SU(3) formalism, as envisioned in [|J. 

In the present paper the low energy spectra and B(E2) transition strengths 
of 21 Ne, 23 Na and 25 Mg are studied using the quasi SU(3) basis and a 
schematic Hamiltonian. The band structure is recovered and the wave func- 
tions along each band are presented. Results are compared both with the 
experimental data and with full shell model calculations ||. Section 2 re- 
views the essentials of the SU(3) model, including a short description of the 
SU(3) basis and the Hamiltonian. In sections 3, 4, and 5 the energy spec- 
tra, band structure and B(E2) transition strengths of 21 Ne, 23 Na and 25 Mg, 
respectively, are presented. Section 6 contains the conclusions. 



2 The SU(3) model 



For a general review of the SU(3) model and the operator expansion in terms 
of SU(3) tensors we refer the reader to references || [14], [16], [l7| , as well as 
to the first article in this series [12]. In what follows we introduce the SU(3) 
basis and the schematic Hamiltonian used in this work. 
The basis states are written as 

[fu]au(K, Liu)}p(\ l^kLlSn, S U }S; JM) (1) 

where is the number of valence protons in the sc?-shell and [/„.] is the 
irrep of the U(2) spin group for protons, which is associated with the spin 
S n = (/* — /^)/2. The SU(3) irrep for protons is (\ w ,fi w ) with a multiplicity 
label cn^ associated with the reduction from U(6). Similar definitions hold 
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(4.0) 28 

(2.1) 1 16 
(0,2)0 10 



(4.1) 1/2 36 

(2.2) 1/2 24 

(3.0) 3/2 18 
(0,3)3/2 18 

(1.1) 1/2 9 



Table 1: Irreps, spins and C2 values for protons and neutrons in 21 Ne. 

for the neutrons, labeled with v. There are other two multiplicity labels: 
p, which counts how many times the total irrep (A, /i) occurs in the direct 
product (A^, /Xt,-) ® (X u ,fi u ) and K, which classifies the different occurrences 
of the orbital angular momentum L in (A,/i). 

The vector states ([J) span the complete shell-model space within only one 
active (harmonic oscillator) shell for each kind of nucleon. As an example 
we take 21 Ne. It has two protons (n n = 2) in the s<i-shell, which can be 
accommodated in three possible irreps: (A^,/!^) = (4,0), (2,1) and (0,2). 
The first and third irreps have spin zero, the second one has spin 1. Each 
one occurs only once (a^ = 1). For three neutrons in the s<i-shell there are 
five irreps: (\ v ,fi v ) = (4,1), (2,2), (3,0), (0,3) and (1,1). Three of them 
have spin 1/2, the other two have spin 3/2. The SU(3) irreps are ordered 
by decreasing values of the expectation value of the second order Casimir 
operator, C2, 



Proton and neutron irreps in 21 Ne are listed in Table 1 with their spin and 
C2 values. 

Calculations performed using the full proton-neutron coupled SU(3) Hilbert 
space in the s<i-shell |Tl|, Q show that the quasi SU(3) [|l| truncation scheme 
is quite efficient. The SU(3) basis is built by taking the direct product of 
the proton and neutron irreps with the largest C2 values and S = and 1 
(for even number of nucleons) or 1/2 and 3/2 (for odd number), and keeping 
from this list only those states with the largest total C2 values [jllj, [12|, [Tj|. 
The proton and neutron representations of each studied nuclei are shown 
in Tables 1, 5, and 8, and the truncated list of their final couplings, which 
describes the Hilbert space used for each nuclei, can be seen in Tables 3, 6, 
and 9. 



((A,/i)|C 2 |(A,/i)> = (A + p + 3) (A + //)-A/i. 



(2) 
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Nucleus 


X 


G-k — G v 


a 


b 


A 


21 Ne 


0.1063 


0.4524 


-0.200 


-0.020 


0.010 


23 Na 


0.0914 


0.4130 





-0.005 





25 Mg 


0.0795 


0.3800 


-0.130 


0.090 


0.024 



Table 2: Hamiltonian parameters in MeV. 



The Hamiltonian used is 



H Hgp -K -\- H S p V ^ X Q ' Q Hp a i r n (3) 

— G v Hp a i r ^ v + a Kj + b J 2 + A sym Ci- 

where H sp%a is the spherical Nilsson Hamiltonian for a = n or v and the 
quadrupole-quadrupole and pairing interaction strengths x> G v and G v has 
been fixed from systematics as in previous work fpTO], |12^ . The parameters 
a, b, and A sym correspond to the three "rotor terms" . They are small, and 
provide some freedom to perform a nuclei by nuclei best fit. The Hamiltonian 
parameters used in this work are listed in Table 2. 
For the spherical Nilsson Hamiltonian 

3 -> 

H sPja = fko {(fj + -) -2k I ■ s- Kftr} (4) 



the Ring and Schuck [19] parametrization k = 0.08 and fi = 0.000 was chosen. 
While both the single particle Hamiltonian and the pairing terms induce a 
mixing of SU(3) irreps, for the values used here the single particle terms are 
by far and away the driving force behind the mixing. 

Hamiltonian (0) proved to be very powerful in the description of the 



normal parity bands in heavy deformed A-odd nuclei p0| , where protons 
and neutrons occupy different major shells. In the description of light nuclei, 
the same Hamiltonian is missing the proton-neutron pairing term, and for 
this reason it is not isoscalar. Given that the main goal of the present work is 
to assert the validity of the quasi SU(3) truncation scheme, we have kept this 
Hamiltonian, using its rotor terms to partially compensate for the missing 
terms. Further applications of this model to light and medium mass nuclei 
must use an isospin invariant Hamiltonian if the model is intended to display 
its full predictive power. 
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(\ n , // 7r )jS' 7r 


(\, Hu)Su 


(A, fi)S total 


(4,0)0 


(4,1)1/2 


(8,1)1/2 


(6,2)1/2 (7,0)1/2 (4,3)1/2 


(4,0)0 


(2,2)1/2 


(6,2)1/2 


(4,3)1/2 


(2,1)1 


(4,1)1/2 


(6,2)1/2,3/2 


(7,0)1/2,3/2 (4,3)1/2,3/2 


(4,0)0 


(3,0)3/2 


(7,0)3/2 




(4,0)0 


(0,3)3/2 


(4,3)3/2 




(2,1)1 


(2,2)1/2 


(4,3)1/2,3/2 




(0,2)0 


(4,1)1/2 


(4,3)1/2 





Table 3: The 17 irreps used in description of Ne. In some couplings the 
total irrep can have spins of 1/2 and/or 3/2. 



In the following sections a detailed analysis of the nuclei 21 Ne, 23 Na and 
25 Mg is presented. The energy spectra, band structure and B(E2) transition 
strengths are shown in each case. 

3 21 Ne 

21 Ne has two protons and three neutrons occupying the sd valence shell. As 
discussed in the previous section, they can be in any of the three SU(3) 
proton irreps and five SU(3) neutron irreps listed in Table 1, which include 
proton states with spin and 1, and neutron states with spin 1/2 and 3/2. 
The truncated proton-neutron coupled basis is listed in Table 3. It includes 
the 17 coupled irreps with the largest C*2 values. 

The 21 Ne energy spectra obtained with this basis using the Hamiltonian 
parameters shown in the first row of Table 2 are presented in the right hand 
side column of Fig. p]. They are compared with the experimental results 
shown in the second column, and with the theoretical results obtained 
with full shell model calculations using the unified s - d shell Hamiltonian 
[ EH] ] (third column). We also show the band structure in Fig. Ej and their 
wave function components in Fig. |3|. 

It can be seen from Fig. p] that the low-lying energy spectra of 21 Ne is well 
described using the truncated quasi-SU(3) space, which includes only the 17 
coupled irreps with the largest C2 values, out of the 137 coupled irreps (some 
of them with p multiplicity larger than 1) which form the full SU(3) basis. 

Despite the very schematic form of Hamiltonian (Q) and the missing 
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E[MeV] 
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6 



5 - 



2 - 



- 



11/2+ 



9/2 5 
9/2' 



7/2+ 



4 1/2+: 
3/2+" 

+ 5/2+: 



5/2 

1/2+ 9 / 2 



7/2+- 
+- 



11/2 



13/2+ 



7/2+ 



H/2+ 

5/ 2 + n /o + 



9/2; 



5/2+ 



7/2+ 



3/2+ 
9/2+ 



1/2+ 
7/2+- 



5/2+ 
3/2+ 



7/2+ 



11/2 



+ 



13/2+ 



1/2+ , 
7/2+ 9 / 2 



7/2+ 



9/2+ 9 / 2 + 
9/2+ 3 
5/2+ 



7/2+ 
3/2+ 



11/2? 



+ 



SU(3) Exp. 



_l/2+ 
-7/2+ 

3/2+ 
-5/2+ 
11/2+ 



SM 



Figure 1: Energy spectrum for 21 Ne. The first column displays the present 
results, the second shows experimental data, and the third column display 
full shell-model results [pT| with a unified sd Hamiltonian [Hi . 
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proton-neutron terms mentioned above, the rotor terms allow a fine tun- 
ing of the spectra. The moment of inertia is slightly increased by selecting 
b = -0.02. It corrects the quadrupole moment of inertia I, where 1/(21) = 
(3/2)x « 0.16. Some excited states are moved to higher energies by using a 
= -0.200. Using A sym = 0.010 for the symmetry term enhances the contri- 
butions of the irreps with A and n even relative to the others, because they 
belong to different symmetry types of the intrinsic Vierergruppe D 2 p2 |. 



It is worth mentioning that, as it was the case in even-even nuclei [12 



the main features of the spectra can be obtained by setting all the rotor 
parameters to zero. On the other hand, if they are values that are too large 
(in absolute value), the whole band structure is destroyed, and with it the 
agreement with the observed B(E2) values (see below). 

The present model predicts the energies of the ground state band, and 
the states 1/2 1; 5/2 2 , 5/2 3 , ll/2 3 , 9/2 3i6 , 7/2 5 , and l/2 2 pretty close to their 
experimental counterparts. The ground state band staggering effect forces 
the states to be clustered by pairs: (3/2, 5/2), (7/2, 9/2), (11/2, 13/2) (see 
also Fig. 2). It is slightly exaggerated in the last case. The levels 3/2 2 , 
7/22, 9/22, 7/24, and ll/2 2 are predicted at energies around 2 MeV lower 
than observed. This is a clear limitation of the model, which can be related 
to the limited and schematic nature of Hamiltonian. Further investigations 
including proton- neutron pairing would clarify this point. In general, the 
results reported here represent a clear improvement from previous SU(3) 
based descriptions . 

Effective charges e n e ff = 1.56e, e u e ff = 0.56 were used in the evaluation 
of B(E2) transition strengths for the three nuclei. B(E2) values calculated 
with the present model, those obtained by the shell model |ZT] and exper- 
imental ones for 21 Ne are shown in Table 4. The agreement with both the 
experimental and shell model values is in general very good. All but the last 
transition connect states belonging to the ground state band. The factor 5 
deviation found in the 5/2f — > 7 /2f transition seem to reflect effects of the 
truncation of the Hilbert space. 

The 21 Ne band structure of this nuclei is presented in Fig. Results for 
the the ground state and K = 1/2 bands are very similar to those found in 
previous studies pi . The staggering in the ground state band noted above 



is clearly seen. 

The SU(3) wave function components of each band are shown in Fig. [3[ 
The percentage different SU(3) irreps contribute to each state is plotted as 
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B(E2) T 


Expt. 


SU(3) 


SM \21\ 


3/2+ - 


-5/2+ 


1.239 ±0.155 


1.215 


1.126 


5/2+- 


-7/2+ 


0.505 ±0.183 


0.710 


0.738 


7/2+ - 


-9/2+ 


0.387 ±0.215 


0.348 


0.391 


9/2+- 


- 11/2+ 


0.248 ±0.165 


0.239 


0.247 


11/2+ - 


- 13/2+ 




0.187 


0.136 


13/2+ - 


- 15/2+ 




0.129 


0.078 


15/2+ - 


- 17/2+ 




0.077 


0.077 


17/2+ - 


- 19/2+ 




0.040 


0.034 


3/2+- 


-7/2+ 


0.675 ±0.151 


0.688 


0.633 


5/2+- 


-9/2+ 


0.906 ±0.097 


0.771 


0.728 


7/2+- 


- 11/2+ 




0.618 


0.707 


9/2+ - 


+ 13/2+ 




0.540 


0.564 


11/2+ - 


- 15/2+ 




0.413 


0.422 


13/2+ - 


- 17/2+ 




0.264 


0.318 


15/2+ - 


- 19/2+ 




0.182 


0.215 


5/2+- 


- 7/2+ 


0.011 


0.002 





Table 4: B(E2) transition strengths for 21 Ne in [e 2 b 2 x 10 2 }. 
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E[MeV] 



13/2+ 



15/2+ 



13/2+ - 
11/2+ - 

9/2+ - 
7/2+ ■ 

5/2+ ■ 
3/2+ 



ll/2 3 + 



9/2+ 



7/2+ 



5/2+- 



3/2+- 



1/2+ 

A = 1/2: 



13/23 



11/2+ 



9/2+ 
7/2+ 



A = 7/2] 



g.s.b. 



Figure 2: Band structure in 21 Ne. 
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1/2+ 3/2+ 5/2+ 7/2+ 9/2+ 11/2+ 13/2+ 15/2+ 



Figure 3: Wave function components of states belonging to a) the ground 
state band, b) the K = 1/2 band and c) the K = 7/2 band in 21 Ne. 
The percentage each irrep contributes is shown as function of the angu- 
lar momentum. The convention used is O for (8,l)l/2[(4,0)0 <g> (4,1)1/2], 
+ for (4,3)3/2[(2,l)l <g> (2,2)1/2], □ for (7,0)3/2[(4,0)0 ® (3,0)3/2], x for 

(6.2) l/2[(4,0)0 <g> (2,2)1/2], A for (6,2)3/2[(2,l)l ® (4,1)1/2] and * for 

(4.3) 1/2[(2,1)1 ® (2,2)1/2]. 
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a function of the angular momentum of the members of the band. Different 
irreps are represented by different lines and symbols (listed in the figure 
caption). All irreps which contributes more than 2% are plotted. The bands 
are recognized by their large B(E2) transition strengths. The slow (adiabatic) 
change in their SU(3) content helps to confirm the band assignment. 

Insert a) of Fig. |3] shows the dominance of the irrep (8,1) with spin 1/2 
in the ground state band up to J=9/2, in agreement with previous studies 
[|23| . [24J. At J = 11/2 there is a clear change in the wave function, which for 
J=13/2, 15/2 is dominated by the (6,2) irrep with spin 3/2. Irreps with spin 
S=3/2 were not included in [p3[| . In the present contribution we are showing 
for the first time its relevance in the low energy spectra. This is the main 
new feature in the present quasi SU(3) truncation scheme. In El it was 



found that even in presence of strong SU(3) breaking interactions, like the 
one proposed by Preedom and Wildenthal, there is a clear dominance of the 
(8,1)1/2 irrep for the ground state, in complete agreement with the present 
results. In the other two bands, insert b) for K = 1/2 and c) for K = 7/2 
in Fig. |3], similar features can be found. While the K = 1/2 band exhibits 
strong mixing between the (8,1) 1/2 and (6,2) 1/2 irreps, the K = 7/2 band 
is dominated by the (6,2) 3/2 band, underlining once again the crucial role 
played by spin 3/2 irreps. 



4 23 Na 

23 Na has 11 protons and 12 neutrons. The valence space contains 3 protons 
and 4 neutrons in the s<i-shell, allowing for the five proton SU(3) irreps and 
ten neutron SU(3) irreps listed in Table 5. From them, only three proton 
irreps and five neutron irreps contribute to the truncated basis, which include 
the 20 coupled irreps with spin S= 1/2, 3/2 listed in Table 6. The complete 
space contains 670 coupled irreps plus their external multiplicities. Notice 
that most irreps can have both spins, due to the 1 $5 1/2 coupling. There 
are also two (7,2) 1/2 irreps, one coming from the (4,1)1/2 ®(3,1)0 coupling, 
and the other from the (4,1)1/2 ® (3,1)1 coupling. 

The right hand side of Fig. f| shows the results obtained for the low 
energy spectra of 23 Na, calculated using the Hamiltonian parameters listed 
in the second row of Table 2, and the Hilbert space described above. This is 
compared with the experimental data [50], presented in the second column, 
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(A-n-, /iyi^SV 


c 2 


(Aj,, Hu)S u 


c 2 


(4,1)1/2 


36 


(4,2)0 


46 


(2,2)1/2 


24 


(5,0)1 


40 


(3,0)3/2 


18 


(2,3)1 


34 


(0,3)3/2 


18 


(0,4)0 


28 


(1,1)1/2 


9 


(3,1)0,1 


25 






(1,2)1,2 


16 






(2,0)0 


10 






(0,1)1 


4 



Table 5: Irreps, spins and C 2 values for protons and neutrons in Na. 



(A.,/,.) 




(A, jj) total 


(4,1)1/2 


(4,2)0 


(8,3)1/2 


(9,1)1/2 (6,4)1/2 (7,2)1/2 


(4,1)1/2 


(5,0)1 


(9,1)1/2,3/2 


(7,2)1/2,3/2 


(4,1)1/2 


(2,3)1 


(6,4)1/2,3/2 


(7,2)1/2,3/2 


(2,2)1/2 


(4,2)0 


(6,4)1/2 


(7,2)1/2 


(4,1)1/2 


(3,1)0,1 


(7,2)1/2,3/2 




(2,2)1/2 


(5,0)1 


(7,2)1/2,3/2 




(3,0)3/2 


(4,2)0 


(7,2)3/2 





Table 6: The 20 irreps used in the description of 23 Na. 
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Figure 4: Energy spectrum of 23 Na, with the same convention of Fig. |T]. 



The SM data are from Glasgow group [25 



14 



E\Me\ 

16 
14 

12 V 17/2+ 15/# 



19/2+_ 

19/2+ 17/2+ 



3/2+ 



19/2+_ 17 ^ 2 + 
17/2+ 15/2+ 



13/2+ 



10 I- , + 13/2+ 15/2^ 13/2+ 

fi L ,, /9+ ll/M ld lf?_ ' ^ = 9/2+ 

13 / 2 i , x 9/2+ 7/2+ 

, . 9/2 3 + H/2+ 5/2+ 

11 ft 7/2 3 + 9/2+ I 1 , 2 * 3/2 3 + 

9 ft 5/2+ 7/2+ ^ 1/2+ 

7/2+ 3/2+^ - j^ = 5/ 2 
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g.s.b. J< " = V 2 

Figure 5: Band structure in 23 Na. 



and with full shell model calculations performed by the Glasgow group [25 



(third column). The ground state, as well as some excited states are well 
described using the quasi SU(3) truncation scheme. On the other hand, as 
can be seen in Fig. |], the model fails to reproduce the energy of the states 
7/2 2 , 5/2 3 , 3/2 3 , 3/2 4 , and others, thereby exhibiting the limits of the model. 
The energy of these states is in general underestimated, probably due to the 
limited Hamiltonian used. 

The band structure, reconstructed from the largest B(E2) calculated val- 
ues, is shown in Fig. [5]. Comparing Figs. |4] and |5] it is noticeable that the 
ground state band is well described, while the excited K=l/2 and K=7/2 
bands are displaced to lower energies than in the observed spectrum. Ad- 
justing the Hamiltonian parameters did not prove to be useful for correcting 
this feature, pointing again to the need of including another terms in the 
Hamiltonian. 

The negative parameter b = -0.005 provides a small correction to the 
moment of inertia, which is similar to what was found for 21 Ne. The other 
two rotor parameters were found to be useless in this nuclei. Both tend to 
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1/2+ 3/2+ 5/2+ 7/2+ 9/2+ 11/2+ 13/2+ 15/2+ 17/2+ 19/2+ 



Figure 6: Wave function components of states belonging to a) the ground 
state band, b) the K = 1/2 band, c) the K = 7/2 band and d) the K = 5/2 
band in 23 Na. The percentage associated with each irrep is shown as function 
of the angular momentum. The convention used is □ for (8,3)l/2[(4,l)l/2 ® 
(4,2)0], + for (9,l)3/2[(4,l)l/2 ® (5,0)1], O for (9,l)l/2[(4,l)l/2 ® (4,2)0], 
x for (7,2)3/2[(4,l)l/2 <g> (3,1)1], A for (6,4)l/2[(2,2)l/2 <g> (4,2)0] and * for 
(6,4)3/2[(4,l)l/2® (2,3)1]. 
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B(E2) t 


Expt. 


SU(3) 


3/27 - 


-> 5/27 


1.573 ±0.233 


1.677 


5/27 - 


-> 7/27 


0.777 ±0.207 


1.123 


5/27 - 


^9/2+ 


1.359 ±0.194 


1.031 


3/27 - 


-> 7/27 


0.987 ±0.085 


1.015 


1/2+ - 


-5/2+ 


0.124 ±0.023 


0.002 


1/2+ - 


- 5/2+ 


2.681 ±0.583 


2.103 


7/2+- 


-9/2+ 




0.447 


7/2+ - 


.11/2+ 




1.216 


9/2+ - 


- 11/2+ 




0.512 


9/2+ - 


+ 13/2+ 




0.888 


11/2+ - 


- 13/2+ 




0.157 


11/2+ - 


- 15/2+ 




1.055 


13/2+ - 


- 15/2+ 




0.322 


13/2+ - 


- 17/2+ 




0.671 


15/2+ - 


- 19/2+ 




0.800 



Table 7: B(E2) transitions for 23 Na in [e 2 b 2 x 10" 



wash out the staggering in the ground state band, and also destroy the band 
structure. For these reasons they were fixed at zero. 

It is remarkable that many states are well described, and that in most 
cases the small deviations with respect to the experimental values show the 
same trend as the shell model results. Some small energy displacements 
from the observed data in the ground state band are similar to those found 
in previous SU(3) studies |26| . 

A study was also performed using the same Hamiltonian but increasing 
the basis for 20 to 302 irreps, around 40% of the full space. The energy 
spectrum showed almost no changes. This is a very strong argument in favor 
of the quasi SU(3) truncation scheme. The most sizable change happened 
for the 5/2i state, whose energy was reduced by 200 keV, and thereby ap- 
proaching the experimental value. 

Fig. ^| shows the wave function decomposition of the states belonging 
to a) the ground state band, b) the K = 1/2 band, c) the K = 7/2 band 
and d) the K = 5/2 band in 23 Na. For each band the percentage associated 
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with different SU(3) irreps is shown as function of the angular momentum. 
The ground state band is dominated by the irrep (8,3) 1/2 up to J = 13/2. 
For larger angular momentum the irrep (9,1) 3/2 has the largest component. 
The K = 1/2 band has a similar structure, but for large angular momentum 
there is competition between the (9,1) 3/2 and (6,4) 3/2 irreps to replace the 
(8,3) 1/2, resulting in a strong mixing. The K = 7/2 band is dominated by 
the (8,3) 1/2 irrep, while the K = 5/2 band is by large made of the (9,1) 3/2 
irrep. 

It must be stressed again that the inclusion of the (9,1) 3/2 and other spin 
3/2 irreps is one of the most relevant features of the quasi SU(3) truncation 
scheme, which accounts for most of its present success, as compared with 
previous SU(3) studies . 

B(E2) transition strengths are presented in Table 7. They are compared 
with the experimental values, available for some transitions between states 
belonging to the ground state band and to the excited K = 1/2 band. The 
agreement is in general very good. The calculated B(E2;l/2i — ► 5/2i) 
value is two orders of magnitud smaller than the experimental one, while 
B(E2;l/2i ^ 5/2 2 ) is larger and well reproduced. 



5 " 5 Mg 

25 Mg has 4 protons and 5 neutrons in the sd-valence shell. The ten proton 
and 12 neutron irreps in which they can be accommodated are listed in Table 
8. From the 1821 proton-neutron coupled SU(3) irreps we selected the 16 
with the largest Ci values to build the truncated Hilbert space. They are 
shown in Table 8, together with their spin. As in the previous cases, some 
coupled irreps can have both spin 1/2 and 3/2. Only the three proton SU(3) 
irreps and the four neutron SU(3) irreps with the largest C*2 values were 
included in the truncated basis. 

The low energy spectra of 25 Mg is depicted in Fig. |7|. The first column 
shows the predicted levels, in the second one the experimental values and 
in third one those found in full shell model calculations. Its band structure 
is shown in Fig. || and the percentage each irrep contributes to the wave 
functions of the ground state, K = 1/2, l/2 2 , 9/2, 13/2 and 11/2 bands is 
presented in Fig. |9|. The Hamiltonian parameters used in this calculation 
are shown in the last row of Table 2. 
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c 2 


(Aj/, /i u )S u 


c 2 


(4,2)0 


46 


(5,1)1/2 


49 


(5,0)1 


40 


(2,4)1/2 


46 


(2,3)1 


34 


(3,2)1/2,3/2 


34 


(0,4)0 


28 


(4,0)1/2 


28 


(3,1)0,1 


25 


(1,3)1/2,3/2 


25 


(1,2)1,2 


16 


(2,1)1/2,3/2 


16 


(2,0)0 


10 


(0,2)1/2,5/2 


10 


(0,1)1 


4 


(1,0)3/2 


4 



Table 8: Irreps, spins and C 2 values for protons and neutrons in Mg. 





(A„,/v) 


(A, 


H) total 




(4,2)0 


(5,1)1/2 


(9,3)1/2 


(10,1)1/2 


(7,4)1/2 


(4,2)0 


(2,4)1/2 


(6,6)1/2 


(7,4)1/2 


(4,7)1/2 


(5,0)1 


(5,1)1/2 


(10,1)1/2,3/2 






(4,2)0 


(3,2)1/2,3/2 


(7,4)1/2,3/2 






(5,0)1 


(2,4)1/2 


(7,4)1/2,3/2 






(2,3)1 


(5,1)1/2 


(7,4)1/2,3/2 






(2,3)1 


(2,4)1/2 


(4,7)1/2,3/2 







Table 9: The 16 irreps used in description of 25 Mg. 
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E[MeV] 



6 - 



5 - 



2 - 



1 - 



- 



11/2+- 

1/2+ - 
9/2 3 + - 



13/2+ 
11/2+ 



9/2+ 



7/2+ • 



5/2+ 
7/2+ = 

1/2+ _ 

9/2+ 



3/2+ 

1/2+ 
5/2+ • 



11/2 



+ 



11/2? 



+ 



1/2+ 



SU(3) Exp 



l/2 3 + 



11/2+ ■ 
= 13/2+ 



11/2+ J 

9/2 3 + 



9/2^ 



+ 



9/2? 



7/2+ 
+ 



- l/2 2 

- 5/2+ 

"7/2+ 

- 3/2+ 

- 1/2+ 



5/2+ -I 



SM 



Figure 7: Energy spectra of 25 Mg. The first column shows the SU(3) results, 
the second the experimental levels, and the third the results of full shell model 
calculations |27| . 
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E[MeV 



16 




14 
12 


15/2+ 


10 


" 13/2+ 


8 
6 


. 11/2+ 


4 


- 9/2+ 


2 


" 7/2+ 





- 5/2+ 



13/2+ 



17/2+ 



19/2+ 



H/2+ 

_15/2 fl + 
13/2+ 



g.s.b. 



H/2+ 15/2+ 17/2+ 11/9+ 

n /2 3 + 13 / 2 | _ _^± /2 < 

15/2+ 

H/2+ 1 3/2+ A' = 11/2 
9/2+ 




A' = 13/2 



K = 9/2 



Figure 8: Band structure of 25 Mg. 

The energy levels shown in Fig. |7| exhibit a complicated pattern. To 
identify the band structure depicted in Fig. ||, both the B(E2) values and 
the wave function decomposition were employed. This band structure is 
fully consistent with the one reported in |28|]. While good agreement has 
been found for a number of states, others like 9/2 1; l/2 2 , 11/2 1; and 13/2! 
show a deviation from the experimental and shell model values, reflecting a 
limitation of the present model, probably due both to the Hamiltonian used 
and the truncation of the Hilbert space. 

The description of the 25 Mg energy levels and its wave functions in terms 
of SU(3) irreps reported here is very similar to the one described in the orig- 
inal work of Draayer p9[j , despite the schematic interaction and the smaller 
basis used in the present work. It strongly supports the reliability of the 
quasi SU(3) truncation scheme. 

The general structure of the energy spectra in Fig. [7| reproduces some of 
the observed levels. The limitations of the schematic Hamiltonian used can 



be gauged by the "band shifts" of three excited bands: those with K = 
1/22, 9/2 and 13/2. The effect of other terms in the Hamiltonian, like proton 
- neutron pairing, in shifting these bands in the right direction will be the 
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Figure 9: Wave function components of states belonging to a) the ground 
state band, b) the K = 1/2 band, c) the K = l/2 2 band, d) the K = 9/2 band, 
e) the K = 13/2 band and f) the K = 11/2 band in 25 Mg. The percentage 
each irrep contributes is shown as function of the angular momentum. The 
convention used is O for (6,6)l/2[(4,2)0 <g> (2,4)1/2], + for (4,7)3/2[(2,3)l 
<g> (2,4)1/2], x for (9,3)l/2[(4,2)0 ® (5,1)1/2], A for (10,1)3/2[(5,0)1 <g> 
(5,1)1/2], * for (10,l)l/2[(4,2)0 <g> (5,1)1] and □ for (7,4)3/2 { [(2,3)1 ® 
(5,1)1/2], [(5,0)1 ® (2,4)1/2], [(4,2)0 ® (3,2)3/2] }. 
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B(E2) t 




Exp. 




SU(3) 




21 
i 


1/27" - 


5/27 





024 


± 


001 


0.398 





530 


3/27 - 


-> 5/27 




n 


033 




0.276 





180 


1/27 - 


-> 3/27 





868 


± 


434 


1.249 


1 


440 


5/27 - 


-> 7/27 


1 


621 


± 


289 


1.968 


1 


813 


3/27" - 


-> 5/2+ 





202 


± 


098 


0.835 





510 


1/27" - 


-> 5 /2+ 


2 


345 

kj t: :J 


± 1 


042 


2.323 


2 


310 


1/2+ - 


^5/27 





156 


± 


074 


0.020 





044 


3/2+ - 


^7/2+ 


2 


084 


± 


261 


1.488 


1 


540 


5/2+ - 


-7/2+ 







010 




0.001 





001 


3/2+ - 


-5/2+ 







074 




0.002 





004 


7/2+ - 


-9/2+ 





814 


± 


271 


0.081*) 


1 


362 


5/2+ - 


-9/2+ 





550 


± 


043 


0.432 





400 


7/2+- 


-9/2+ 





054 


± 


038 


1.420*) 






5/2+- 


- 9/2+ 





097 


± 


006 


0.545 






5/2+- 


- 9/2+ 


2 


171 


± 


362 


2.066 






7/2+- 


* 11/2+ 





847 


± 


261 


0.001 






7/2+- 


♦ 11/2+ 





332 


± 


078 


0.372 






9/2+- 


♦ 13/2+ 





128 


± 


061 


0.279 







Table 10: B(E2) transitions for 25 Mg in [e 2 b 2 x 10~ 2 ]. The second column 
shows the experimental values, the third column the present theoretical re- 
sults, and fourth those obtained previously with the SU(3) model p9| . Stars 
denote transitions with unclear identification in the theory. 



subject of future research. 

In view of the limitations of Hamiltonian (^), a very specific selection of 
rotor terms was needed to get a 5/2 ground state. The parameter set shown 
in the last line of Table 2 reflects a subtle balance between the J 2 , K 2 and 
symmetry terms. The negative a value pushes bands with larger K down 
in energy, while the positive b reduces the moment of inertia, increasing the 
level separation inside each band, and moving the band heads with larger J 
to higher energies. Setting A sym = 0.024 helps to make the J = 5/2i state 
the ground state. 

Table 10 lists the B(E2) transitions for 25 Mg, calculated with the same 
effective charges e n e ff = 1.56e, e u e ff = 0.56e as used in the previous cases. 
The second column shows the experimental values, the third column the 
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present theoretical results, and fourth those obtained previously using the 
SU(3) model with a renormalized Kuo-Brown interaction ||29|| . The agree- 
ment is in general quite good. 

Stars were used to mark the calculated transition strengths B(E2; 7/2\ — > 
9/2i) = 0.081 ([e 2 6 2 x 10~ 2 ]) and B(E2; 7/2 1 -> 9/2 2 ) = 1.420. Due a shift 
of the K=9/2 band, the 9/2i and 9/2 2 are inverted in our level scheme. A 
more natural assignment would be B(E2; 7/2 1 — > 9/2i) = 1.420 and B(E2; 
7/2i — > 9/2 2 ) = 0.081, which are close to the experimental results. 

Calculated results for the B(E2) transition strengths are not very good. 
Only a few are well described: 1/2^ -> 3/2f, 5/2^ -> 7/2|, 1/2| -> 5/2^, 
5/2f -> 9/2f , 5/2^ -> 9/2^, 7/2f -> ll/2f , and 9/2f -> 13/2f . For the 
other reported transitions the difference is one or more orders of magnitude, 
most being underestimated. These wrong predictions show the limitations of 
the model. It is interesting to see that the present values and those reported 



in Ref. |29| are in general very close. 

As mentioned above, by using the B(E2) transition strengths and the 
form of the wave function it was possible to identify the six rotational bands 
shown in Fig. || in agreement with |28j. The SU(3) content of these bands 
is shown in Fig. |9] as a function of the angular momentum of the states 
belonging to the different bands. Beyond the regular structure found in 
most of the bands, the most remarkable feature is that of first two bands 
which have band heads that are dominated by the (6,6) irrep. The "leading" 
irrep (9,3)1/2 is dominant in the second excited K=l/2 2 band. This is a 
remarkable result: the quadrupole - quadrupole interaction builds the ground 
state band mostly from the leading irrep. As was pointed out in [p9|] , it also 
implies a coexistence of prolate and triaxial shapes, associated with the (9,3) 
and (6,6) irreps, respectively. It also shows that, while developing clear 
rotational bands, the 25 Mg ground state is mostly triaxial. The yrast band, 
however, includes many different band heads, which are either triaxial or 
prolate. 

At variance from what was found in the lighter nuclei discussed above, 
the spin 3/2 irreps play only a marginal role in 25 Mg, except for the states 
with the largest angular momentum in the K = 9/2 band, insert d), and the 
K= 11/2 band, insert f), which is dominated by the (7,4) 3/2 irrep. 
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6 Discussion and Outlook 



The quasi SU(3) truncation scheme was used in conjunction with Hamilto- 
nian (||]) to describe the energy spectra and B(E2) transition strengths in 
21 Ne, 23 Na and 25 Mg. Comparison were made with experimental data, with 
shell model calculations and with previous SU(3) studies. The agreement 
was in general good, exhibiting the success of the model. 

The truncation recipe is quite simple. First select three or four proton 
and neutron SU(3) irreps with largest C2 values, which in general will have 
spin or 1 for an even number of nucleons and spin 1/2 and 3/2 for odd 
an number of nucleons. Build from these the coupled proton-neutron SU(3) 
irreps, and select again only a small number with the largest total C2 values. 
Usually 20 irreps are enough. Use this basis to diagonalize the Hamiltonian. 
While the effect of the truncation in 21 Ne was small, for 25 Mg it implied a 
two order of magnitude reduction in the basis size. The composition of the 
calculated wave functions was found to be very similar to the ones reported 
in previous SU(3) studies with far larger basis where the Kuo interaction was 
used fX$. 

The Hilbert space built in this way is rich enough to include many excited 
rotational bands which have a clear counterpart both in the experimental 
data and in full shell model calculations. The fact that most of these bands 
have one fairly dominant SU(3) irrep underscores the strength of the model. 
This feature leads to a simple picture in terms of the SU(3) decomposition of 
the wave functions, and exhibits the crossing and mixings which occur within 
a band as a function of the angular momentum of the different band members. 
B(E2) transition strengths were also found to be in close correspondence with 
the experimental data, and this allowed for a clear identification of band 
members. The SU(3) content of the ground state band also showed that, 
while 21 Ne and 23 Na are definitively prolate, 25 Mg is mostly triaxial. 

There were model limitations uncovered in the present study. The most 
stricking are the band shifts, i.e. the fact that some excited bands are pre- 
dicted at energies 1 to 2 MeV lower than they appear in the experiment. 
Given that the band structure and the B(E2) values were in general well 
depicted, the band shifts seem to reflect on a limitation of the Hamiltonian 
and not of the Hilbert space, since the same feature was found in full shell 
model calculations. 

It should be clear that Hamiltonian (H) is not missing important two-body 
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terms. In addition to realistic single particle energies, it includes quadrupole- 
quadrupole and like particle pairing terms with fixed interaction strengths 
taken from systematics. It has also three rotor like terms which allow for 
a fine tuning of the spectra. In the case of 25 Mg, these terms played an 
important role in the reproducing of the energy spectra, especially in pushing 
the J = 5/2 band down to become the ground state band of the system. 
In these nuclei with N very close to Z, proton-neutron pairing, both in 



the T=0 and T=l channels, plays a very important role [30|. Its inclusion 



in Hamiltonian @ would make it isoscalar. It could be relevant not only 
for improving the predicted energy spectra, but also in the description of Ml 
excitations, which are known to be very challenging for any theoretical model 
It would also allow the model to be tested in the fp-shell where B(M1) 
transition strengths have been recently measured [32]. Future research on 
these subjects is desirable. 
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